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EFFECT OF THE ACCELERATION OF ELONGATED BODIES CF RE70LUII0N 
UPON THE RESISTANCE IN A COMPRESSIBLE GAS* 

By F. I. Frankl 


The problem, of the motion of an elongated, "body of re-volution in 
an incompressible fluid, may, sis is known, "be solved, approx i mately with 
the aid. of the distribution of sources along the axis of the "body. In 
determining the velocity field, the question of whether the tody moves 
uniformly or with an acceleration is no factor in the problem. The 
presence of acceleration must he taken into account in determining the 
pressures acting on the body. The resistance of the body arising from 
the accelerated motion may be computed either directly on the basiB of 
these pressures or with the aid of the so— cal led associated masses 
(inertia coefficients). A different condition holds in the case of the 
motion of bodies in a compressible gas. In this case the finite 
velocity of sound must be taken into account. If it is assumed that 
the body produces in the flow only small disturbances, the velocity 
potential cp satisfies the wave equation: 

a3&p _ (l) 

St 2 

where a is the velocity of sound. The method of sources still remains 
applicable but in computing the effect of the sources it is necessary 
to make use of the retarded potential. 

We Introduce two systems of coordinates, one a fixed system xp, 
the other moving with the body where x and I are, respectively, 

the coord inat es along the axis and r is the distance from the axis. 

The coordinate I is computed from the nose of the body. 

Let r = r(i) be the equation of the body of revolution (0 < | < 2) . 
The position of the body is characterized by the inequalities 
— f(t) < x < 2 — f(t), where f(t) is a given function of the time 

characterizing the motion of the body. The coordinates x and I are 
connected by the equation x = | — f (t) . 


*”0 Yliianii uskoreniya na soprotivlenie pri dvizhenii prodolgovatykh 
tel vrashcheniya v gazakh." Erikladnaya Matematika i Mekhanika, Tol. 10, 
no. 4, 19 ^ 6 , pp. 521-524. 
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Let q = q(x’ ,t’) be the value of the linear densi ty of the source 

at the point x* at the instant of time t* and B = \J (x — x’ )^ + r^ 
he the distance between the source and the point considered. The i 
potential of the disturbance velocities is then 

where the integral is taken between the limits determined by the 
inequalities 


-f 



< x* < l - f 



In order to establish the range of integration, we must consider 
in the plane x*t* a strip — f(t*) < x* < 2 — f (t 1 ) between the 
lines of motion of the front and rear points of the body (fig. l) and 
the part of the semihyperbola in this strip. 

t* = t — S = t — i\Z(x’ — x) 2 + fB 2 2 - ~ 

a a v wy 


The abscissas of the part of the semihyperbola considered form the 
range of integration (fig. l) . As may be seen from the figure, the 
range of integration may consist of one or several segments. The ends 
of the corresponding arcs of the hyperbola may lie either on the line 
of motion of the nose or on the line of motion of the base. In 
particular, in the motion of the body with supersonic velocity it is 
possible that within the strip considered there lies an entire arc of 
the hyperbola, both ends of which lie on the line of motion of the 
nose (as will be the case, for example, for a uniform supersonic motion 
or for a motion approximating this type). 

The function q(x* ,t‘) must be determined from the integral 
equation expressing th® fact that the normal component of the velocity 
of the gas on the surface of the body is equal to the normal component 
of the motion of the body: 


5cp 

5n 



B 4 (*''* 't 



dx* 


q/« f .(t) 

\J 1 + [r* U3 d 


( 2 ) 


where S/Sn denotes differentiation along the outer normal and 
I = x + f(t). In what follows we consider the derivative dr/di as 
small, and shall give an approximate solution of equation (2). For the 
assumption made, Scp/Bn may be replaced by Scp/Sr. We assume for 
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simplicity that the range of integration reduces to the segment (x^x^). 
Equation (2) then assumes the form 


r 



If r/2 and dr/d| are considered small quantities of the first order 
and their higher degrees or derivatives are neglected, an approximate 
solution of equation (3) will be 1 


q(x,t) = 2itf * (t)r(j-) 


w 


¥e shall show that the error in substituting this expression in 
equation (3) is a magnitude of the second order of smallness. Eor this 
Impose it is convenient to write er* for r where r* = r*( |) is a 
variable quantity of the order of unity relative to l and e is a 
small constant number. The degree of e in each expression then shows 
its order of smallness. We consider the order of the terms of the left 
side of equation (3). The terms outside the integral on the left of 
equation (3) containing q as a factor for finite terms (in the general 
case) will be of the second order of smallness. In the same way the 
expression 

P*2 . 

r / oq dx* 

E 2 


<Me 2 r 



(5) 


is of second— order smallnesB. We write the remaining term in the form 




■HChis solution (4) was first used by Karman (reference 2) in the case 
of a uniform motion ( in particular, for obtaining the approximate value of 
the wave resistance) . 
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As a result of the computation of the first integral of this expression 
■we obtain 


r 



1 

er*(t) 


2 + 0(e 2 )] 


( 7 ) 


For the estimates 
and third integrals of 


of the differences in the brackets in the second 
expression (6) we have, respectively 


q.(x*,t) - < 3 .(xj t ) = 0(e 2 )E, qXx’^t) - q(x*,t - E/a) = 0(e 2 )B (8) 


Hence the estimates of integrals are of the form 


0(e 2 )r 



rt0(e 2 ) 

E 2 - 


(9) 


so that these integrals are likewise magnitudes of the second order of 
smallness . Thus 


1 A 

Sr 


f L 


dx* = f*(t) ^ + 0(e 2 ) 


( 10 ) 


as was required to be proven. 

The egression for the velocity potential after substituting in 
equation (4) becomes 

cp(x,t) = fi f* (t -§) 7(|») cbc* (V(|) = (11) 


For computing the field of pressures we have the generalized formula 
of Bernoulli 



( 12 ) 


where p is the pressure at the given point, p the pressure in the 
undisturbed 'region, p the density and, w the modulus of the velocity. 
Neglecting magnitudes of the second order w 2 /2 and taking account of 
the equations 


1 ' = x- + f (t - 1) , 


Ml = f ft - S) 

St \ a / 


(13) 
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we obtain 



In the particular case where the pressure is required, on the surface 
of a projectile, the motion of which approximates to uniform supersonic 
motion, the point Xq lies on the line of motion of the nose of the 

projectile and the last term on the right of equation (lh) is equal 
to zero. We obtain 



where x-j_ and, x^ are determined from the equation 



In the case of uniform motion x = |* — vt + c and equations (15) 
and (l6) become 



The additional pressure produced by the acceleration is therefore given by 
the equation 
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where £*• is determined, from the equation 


X* = £*f _ ft(t)t* + f*(t)t - f(t) 


( 19 ) 


From the structure of equation (18) it is easily seen that the 
relative increase of the pressure arising from the acceleration has, for 
velocities comparable with the velocity of sound, an order of magni- 
tude bl/v 2 where b = f"(t) is the acceleration. Hence for rocket 
projectiles of the usual dimensions, for accelerations not higher 
than 1000 meters per second 2 and velocities comparable with the velocity 
of sound, the added pressures arising from the acceleration are 
negligibly small. 


This is confirmed also by the following illustrative computation. 

The length of the war head portion of the shell is equal to 2 = 0.25 meter 
and the maximum radius r = 0.07 meter. The generator of the head is an 
arc of a parabola which smoothly goes over into the cylindrical part so 
that it is given by the equation r = O.56I — 1.12| 2 (scale in meters). 

The motion is one of uniform acceleration with an acceleration 
of 1000 meters per second 2 so that the line of motion of the 
projectile x = — 5000t 2 (t > 0), x = 0 (t < 0) . 


The pressure distribution over the war head is found at the instant 
of time t = 0.5 second such that the velocity is equal to v = 500 meters 
per second. The air density p = 0.125 kilograms - Iper second 2 . There 
were obtained for the pressure distributions (with account taken of the 
acceleration) and the additional pressures 5p produced by the 
acceleration the following values: 


I = 0.05 0.10 0.15 0.20 i = 0.05 0.10 0.15 0.20 (meters) 

p - p = 8550,7448,5970,3588; 5p = -5.0 -5.5 -11.0 -16.0 (kilograms per meter 2 ) 


The wave resistance and the corresponding added resistance were then 
obtained as 



(P - P) 


cLS = 100, 


r j£sp d£ = 0.13 
d§ 


Q = 2rt 


H 


5Q, = 2n 


0 


(kilograms) 

( 2 0 ) 
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The added pressures Sp arising from the acceleration were computed from 
the appro -si mate equation (l8) which for the example given has the form 



where 


Sx-[_ = xj_ — X 2 * & 



bt Z 1 “ 1 9 

a 


5xg = Xg — Xg* ~ 



( 22 ) 


Translated by S. Eeiss 
National Advisory Committee 
for Aeronautics 


1. Lamb, H. : Hydrodynamics . Cambridge, 6th edition, 1932. 

2. v. Barmin, Th. : Problem of Eesistance in a Compressible Fluid. 1935* 
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Figure 1. 


HADA TM 1230 



